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DECUS Program Library Write-up 


DECUS NO. 8-485 


1. Program Title : Geometric Data Truncation for Fourier 

Transform Programs 

2. Abstract : This program is intended for use as a truncating- 

weighting subroutine in conjunction with a Fourier 
Transform program. The application of a weighting 
function to the data record before the application 
of a Fourier transform program reduces the spread in 
frequencies which results from the transformation 
of a finite record. This procedure is often called 
apodization in that it removes the side lobes in the 
transformed record that results from a rectangular 
data window. 

3. Requirements : 

3.1 Storage : The subroutine occupies 115(8) locations and 
is relocatable at the time of assembly. 

3.2 Hardware : PDP-8 

3.3 Software : The user must have a mainline program for 
gathering and storing in serial order the data before 
the application of this truncation algorithm. Following 
the weighting of the data the transform must be applied. 
The real series must be stored in serial order and 
without gaps. The length must be of the form 2N, just 
as required by the FFT. Additionally, the data must 

be in signed 11 bit, fixed-point format. 

4. Usage : The program is intended to be assembled along with 

the mainline program. It is to be called after the 
data are gathered and before the actual FFT. The weighted 
data are replaced in their original locations with 
the original'Contents being lost. The subroutine is 
called with the number of data points in the AC and the 
address of the first point in CALL +1. Return is to 
CALL +2. 

5. Discussion : When performing the discrete Fourier transform 

on a finite record in the time domain an inherent 
distortion results in the frequency domain. 

The FFT gives the power spectral density of an observed 
signal. Denoting Fourier transforms by ( ) and 

convolutions by (*), 

A(v) = ATtT $ F (v) = FTET 


A(v) = \ A(t)e"^^ 1TVt dt 
J - 00 


where 
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If F(t) is the temporal representation of the observed 
signal, the power spectral density is given by 


F(v) = 

A°° 

\ F(t)e' l27TVt dt 

= 

fTET 


V - oo 




The data window here extends from - 00 to +°° and weights 
each datum equally, i.e., A(t) = 1 for - 00 ^ t^ +°°. 

The multiplication of two functions in the time domain 
is equivalent to convolving the transforms of the two 
functions in the frequency domain. 


A(t) • F(t) = AXt7*FTtT = A(v)*F(v) 



v ' ) F(v)dv' 


For the case of an infinite data window 


00 

A(v) = AltT = C l-e _i27TVt dt = 6 (v - v’) 


Hence the power spectral density 



A 30 


(I 00 

F(v) = 

C F(t)e'^ llrvt dt 

= 

V 6(v - v ! )F(v)dv 1 


J - 00 


J - 00 


is the true power spectrum. 

Infinite records are impossible and the observations 
of the temporal signal must be terminated. This 
limiting of F(t) is truncation. 

The most common form of truncation is the rectangular 
data window. Observations are made for a certain period 
of time and all information before and after this period 
is neglected. This gives rise to an A(t) that is a 
rectangle of unit height. The Fourier transform, 

A(v) = sin uv , is a function with a finite spread and 
irv 

smaller side lobes, (see Fig 1A). The true spectrum is 
then convolved with a function which spreads each 
frequency into several frequencies. 

Several other truncation functions have been used. 

These functions reduce the spread of the central peak 
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and/or remove the side lobes. Two of these are the 
triangular function with A(v) = / sin ttv ) 2 (see Fig. IB). 

V TTV ' 2 

and the Gaussian function with A(v) = e _Trv , (see 
Fig. 1C). The application of these functions involves 
a degree of data handling that may be overly time 
consuming in order to obtain a clean power density 
spectrum. 

This program involves a geometric truncation which 
requires no multiplications, or divisions, and yet 
gives a smooth truncated signal to be transformed. 

In the frequency domain the function to be convolved 
with the true spectral function has no side lobes 
and a minimal width. 

6. Description: The data record is first divided into 16 

blocks. The binary representation of each datum 
within a block is shifted to the right a number of 
times dependent upon the position of the block in the 
data record. The components of the blocks would be 
shifted 7, 6, 5,..., 1, 0, 0, 1, 2,..., 6, 7 places, 
resulting in multiplications of the components by 
1/128, 1/64, 1/32,...1/2, 1, 1, 1/2, 1/4,..., 1/64, 1/128 
respectively. These factors are indicated by the hori¬ 
zontal bars in Fig. 2. 

The truncation function when smoothed to its average 
value is the hyperbolic secant, sech (irt) . The geo¬ 
metrical function has a number of discontinuities but 
at points of discontinuity the Fourier transform 
integral converges to the average 

limit —[A(t * e) + A(t - e)] 
e+0 2 

of the right and left hand limits. 

The Fourier transform of sech (irt) is sech (ttv) , (see 
Fig. ID)- So the function to be convolved with the 
true spectrum has no side lobes and minimal width. 

7. Execution Time : The time required to weight a data record is 

(100 ps) x (data record length). 

For example, a 512(10) point record length can be 
weighted in 0.0512 seconds. 

8. User Modification : The user who intends to transform data 

records of a fixed length may wish to shorten the 
subroutine by calculating and storing some of the 
constants which would otherwise be recalculated and 
stored at each execution time. 
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9. Listing : A listing is below. 

10. Logic Flow Chart : A logic structure flow chart follows 

the listing. 

11. References : 

(T) Bergland, G. D. "A Guided Tour of the Fast Fourier 
Transform", IEEE Spectrum , 6^, 41 (1969). 

(2) Bingham, C., Godfrey and Tukey, "Modern Techniques 
of Power Spectrum Estimation." IEEE Trans. Audio 
Electroacoustics, AU-15, 56 (1967). 
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Figure 1 


















Truncation Function (Geometric 5 sech ttx) 
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Figure 2 
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/GFNFPAJ. I Z FD AFODIZATION FUNCTION 
/F. A. FAPNHAPPT, 3P SFPTFMPFR 197 1 

/FOP DATA OF LENGTH P»N, THAT IS 10PZ, 5 1 P, P5f, FTC 

/PUNNING TINE , 100 MICPOSEC TINES DATA LENGTH 

/DATA RFSTOFFP TO OFIGINAL LOCATION. 

/CALL WITH NUMFF'F OF DATA FOINTS IN AC AND ADDRESS 
/OF FIRST POINT IN CALL+1, PFTUPN IS TO CALL+P 
/FROGRAM OCCUFIFS 115(8) LOCATIONS 

* /< f pp 

4600 0000 APOD, 0 


/If 01 

3313 


DCA 

PTS 

/If OP 

1 fOO 


TAD 

I APOD 

ZfP3 

331 1 


DCA 

ADD! 

A 60 A 

PPOO 


I S7. 

APOD 

4605 

1313 


TAD 

FTS 

/if Of 

7 0 1 P 


FTP 


Af 07 

7 01 P 


R TF 


Af 10 

330P 


DCA 

TEMP? 

Af 1 1 

1 30P 


TAP 

TEMP? 

Af 1 P 

7 0 A 1 


CIA 


Af 13 

330 f 


DCA 

CN TP 3 

Af 1 A 

1313 


TAD 

FTS 

Af 1 5 

1 30 f 


TAD 

CN TP3 

Af 1 f 

1 30 f 


TAD 

CN TP 3 

Af 17 

7 O/'l 


CIA 


Af PO 

3307 


DCA 

CNTPZ 

Af P 1 

1313 


TAD 

FTS 

Z'fPP 

7010 


PAR 


AfP3 

1311 


TAD 

ADD 1 

AfP/i 

1 30? 


TAD 

TFMPP 

Zf? 5 

33 IP 


DCA 

ADD? 

Af P f 

131 A 


TAD 

M7 

Af P7 

3303 


DCA 

CN TRO 

Af 30 

1 30 f 

L. OOF 1, 

TAD 

CN TI 3 

Af 3 1 

3310 


DCA 

CNTF5 

Af 3P 

17 1 1 

SLOOP 1, 

TAD 

I ADD 1 

Af 33 

zp f 5 


JMS 

SFSR 

Af 3A 

37 1 1 


DCA 

I A DD1 

Af 3 5 

P31 1 


I SZ 

A DPI 

Af 3 f 

P307 


I SZ 

CNTPZ 

Af 37 

P30 f 


I SZ 

CN TP 3 

Af AO 

5P3P 


JMP 

SLOOP 1 

Af A 1 

P30 3 


I SZ 

CNTPO 

/if /i? 

5P30 


JMF 

LOOP 1 

Af A3 

7 300 


CL. A 

CLL 

Af A/i 

7001 


I AC 


/if A 5 

7 0 A 1 


CIA 


AfAf 

3303 


DCA 

CNTPO 

Af A7 

1 30 f 

LOOP?, 

TAD 

CN TP 3 

Zf 50 

3310 


DCA 

CN TF 5 

zf 51 

17 IP 

SLOOFP, 

TAD 

I ADD? 

Zf 5? 

A? 65 


JM S 

SF SP 

Zf 53 

37 IP 


DCA 

I ADD? 

zf 5/- 

P31P 


I SZ 

ADD? 


/NUNPEP OF DATA FOINTS (8) 

/ADDRESS OF FIRST DA TIN 
/SF T RETURN TO CALL +P 


/L.FNGTH OF ONE DATA FLOCK 


/MINUS FLOCK LENGTH 


/ACTUAL # OF POINTS TO PF SHIFTED 


/ADDRESS OF PND FLOCK OF PND HALF 
/MINUS 7, MINUS f, MINUS 5, ... 


/GET DATUM 
/SHI FT 

/PFPOSIT SHIFTED DATUM 
/UPDATE ADDPFSS 
/ALL COMPLETE? 

/NO, FLOCK COM FL FTP? 

/NO, GET NEXT DATUM 

/YFS, REDUCE SHIFT OP GO TO PND HALF 

/DO NFXT PLOCK 

/SF CON D HAL F INI TI ALI 7. AT I ON 


/MINUS 1, MINUS P, MINUS 3, ... 


/OFT DATUM 
/SHI FT 

/DEPOSIT SHIFTED DATUM 
/UPDATE: ADDRESS 
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46 5 5 

P307 


I sz 

CNTPA 

4656 

5P60 


JMP 

SKIP 

4657 

5600 


jmp 

I A POP 

4660 

P310 

SKIP, 

I sz 

CN TP R 

46 61 

5P51 


tJM F 

SLOOPP 

466P 

1303 


TAP 

CNTPO 

4663 

7 04 1 


CIA 


4 664 

5 P 47 


JT*P 

L 00 FP 

4665 

0000 

SPSP, 

0 


4 666 

3301 


ncA 

TFMF 1 

4 6 67 

1 303 


TAP 

CMT r -0 

467 0 

3 3 1 5 


pca 

CM T 

467 1 

1 30 1 


TAP 

TEMPI 

4 67 P 

7 100 

L00F3, 

DLL 


467 3 

7 510 


SPA 


467 4 

7 0 P 0 


CM, 


4 67 5 

7 010 


PAP. 


467 6 

P3 1 5 


I SZ 

CN T 

467 7 

5P7P 


p 

L.00P3 

47 00 

5 6 65 


r 

I SPSP 

47 0 1 

0000 

TFKJ 1, 

0 


47 0P 

0000 

TFMFP, 

0 


47 0 3 

0000 

CM TP 0, 

0 


47 04 

0000 

CMTP ! , 

0 


47 0 5 

0000 

CNTFP, 

0 


47 0 6 

0 0.0 0 

CMTR3, 

0 


47 07 

0000 

CM TP A, 

0 


47 1 0 

0 0 0 0 

CM TP R, 

0 


47 11 

0000 

Am, 

0 


47 IP 

0 0 00 

A PUP, 

0 


47 13 

0000 

PTS, 

0 


47 1 4 

777 1 

M7, 

777 ] 

l 

47 15 

0000 

CM T, 

0 


APD1 

47 1 1 




appp 

47 1 P 




PFOD 

4 60 0 





CN T A 7 l 5 

CNTRO A 7 P 3 

CNTF 1 41 Pit 
CNTFP 41 p 5 

CNTR3 47 Ef 

CNTpA /i7 07 

CN TP R 47 1 (’ 

LOOP 1 4 esc* 

L'OCFP A (47 

LC'OP? A ei9 

M7 in 1 4 

RTS 47 1 3 

SKIP 4f(C 

SLOOP 1 4 f 3P 

SLOOPP Af R 1 

SFSF A 6f R 

TFMPl A 7 p l 

TEMPP A7 0 P 


/ALL COM PL F TP ? 

/NO, SKIP 
/YF S, RETURN 
/FLOCK COMFLETE? 

/MO, GFT NEXT DATIM 

/YFS, INITIALIZE FOP NFXT FLOCK 


/SHIFT SUP. (MODI El ED RFC- 08-EM JA- D) 

/STORE NIMFFP 

/(- F T SHI F T COUN I FF 

/STUFF FOP LOCAL USE 
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Yes 


1 
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Yes 


1 
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